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Abstract Measurements of inclusive observables, such as particle multiplicities and momentum spectra, 
have already delivered important information on soft-inclusive ("minimum-bias") physics at the Large 
Hadron Collider. In order to gain a more complete understanding, however, it is necessary to include also 
observables that probe the structure of the studied events. We argue that forward-backward (FB) correlations 
and event-shape observables may be particulary useful first steps in this respect. We study the sensitivity of 
several different types of FB correlations and two event shape variables — transverse thrust and transverse 
thrust minor — to various sources of theoretical uncertainty: multiple parton interactions, parton showers, 
colour (re)connections, and hadronization. The power of each observable to furnish constraints on Monte 
Carlo models is illustrated by including comparisons between several recent, and qualitatively different, 
Pythia 6 tunes, for pp colHsions at ^/s = 900 GeV. 

Key words. Inelastic hadron-hadron collisions - Minimum bias - Forward-backward correlations - Trans- 
verse thrust - Multiple parton interactions - Event generators - PYTHIA tunes 



X 
J3 



PACS. XX.XX.XX No PACS code given 



1 Introduction 



Until recently, the measurements used to constrain physics 
models of high energy particle collisions came primarily 
from experiments done at the previous generations of accel- 
erators, such as the SPS, LEP, and the Tevatron. In particu- 
lar, studies of "minimum-bias" and underlying-event physics 
have been widely used to constrain the poorly known non- 
factorizable and non-perturbative aspects of Monte Carlo 
(MC) event generators. These generators are, in turn, used 
ubiquitously over a continually expanding range of ener- 
gies and intensities, for both high- and low-p_L processes. 
MC "tunes" that rely exclusively on these older data sets 
are, however, becoming outmoded by a new generation of 
high-energy experiments, performed at the Large Hadron 
Collider (LHC). The extrapolation of previous results to the 
higher energies and large acceptances of the LHC exper- 
iments is associated with significant uncertainties, and the 
demands on both theoretical and experimental precision are 
becoming ever more stringent. The importance of reevalu- 
ating the physics models, and of retuning them in situ, has 
therefore been highlighted in several recent studies lHJQ. 



The LHC offers a rich cornucopia of opportunities to 
test and expand the data sets used for MC tuning. Moving 
from the low-p^ results of minimum-bias to the underlying 
event in hard processes, and from central rapidities to ones 
close to the beam axis, we may test the universality and 
applicability of the modeling over a large dynamical range. 

For minimum bias (MB), which is the focus of this pa- 
per, there is no "hard scale", and hence all observables re- 
ceive large non-perturbative corrections. From the point of 
view of a theoretical modeling based on factorization and 
perturbative QCD, it is difficult to say anything meaning- 
ful about this data set, except perhaps for the tiny fraction 
of it that includes hard jets. In Pythia, as in most other 
contemporary MC models, the modeling of soft-inclusive 
physics, is based on the concept of Multiple Parton Interac- 
tions (MPI). Though these corrections go beyond the reach 
of the standard factorization theorems, they are still regarded 
as essentially perturbative in origin; they are dominated by 
(multiple) t-channel gluon exchanges, regulated by a non- 
perturbative dampening at low px, and dressed with non- 
perturbative descriptions of the beam remnants and of the 
hadronization process. We shall not here go into the de- 
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tails of the modeling (for a recent review, see It is 
worth noting, however, that the smooth transition between 
soft and hard scattering processes is the main reason the 
Pythia modeling can be used for both minimum-bias and 
underlying-event physics. As mentioned above, this univer- 
sality is important, since by testing the model in a minimum- 
bias environment here, we may expect to simultaneously 
improve the description also of other physical processes. 

Particle multiplicities and transverse momentum spec- 
tra are generally used to give the first important constraints 
on the overall amount of particle production and on its dis- 
tribution in and rj. Several such studies have already 
been published by the LHC collaborations ll9 UT7l . How- 
ever, as most people familiar with MC tuning will be aware 
of, there are often several qualitatively different ways of 
mixing the same cocktails. (I.e., several tune parameter sets 
may reproduce the same experimental data.) An important 
question therefore concerns the balance between the sev- 
eral different particle production mechanisms that are avail- 
able to an open-minded model builder: initial- and final- 
state radiation, beam remnant breakup, "hard" processes 
vs. additional MPI interactions, final-state interactions, etc. 
Each production mechanism generally has some kinemat- 
ical or dynamical "signature" that can be used to single it 
out, if sufficiently differential information is available. Tests 
using several, mutually complementary, discriminating ob- 
servables are therefore to be regarded as essential to over- 
come model degeneracies. 

In the context of the present paper, it is especially impor- 
tant to note that the collinear enhancements characteristic of 
parton shower activity tend to produce local I short-range 
correlations, i.e., the additional particle production caused 
by the showering dies away rather quickly with distance to 
the originating parton. By contrast, the particle production 
associated with semi-hard MPI (mini-)jets tends to produce 
an enhancement of correlations around A(j) = vr while soft 
particle production from strings stretched between the rem- 
nants should generate long-range correlations in rapidity 
which, to a first approximation, should be homogeneous in 
(p. Thus, the shape of events and the particle -particle corre- 
lations within them, can be used to gain a handle on the rela- 
tive strengths of different particle production mechanisms. 

In section |2l we briefly present the Monte Carlo mod- 
els we shall use, and discuss the various sources of par- 
ticle production relevant to soft-inclusive physics. A mild 
selection bias is then introduced in section [3j to mimic a 
"minimal" minimum-bias selection, and the effect of in- 
cluding additional cuts is illustrated. In section |4l we 
briefly compare our reference models on some of the typ- 
ical minimum-bias plots, such as charged multiplicity, r], 
and p±^ distributions. In section |5l we then turn to a more 
detailed study of forward-backward correlations, including 
ones with an explicit cj) dependence, and in section|6]we dis- 
cuss the transverse thrust and the transverse minor. We give 
concluding remarks in section |7] 



2 Monte Carlo Models and Parameters 

For details on the modeling of hadron collisions incorpo- 
rated in general-purpose event generators, we refer the reader 
to the recent review JBl. Here, we use the Pythia 6 gen- 
erator throughout and focus on a few main aspects of the 
modeling that it will be useful for the reader to be aware of. 

Firstly, the modeling of the diffractive contribution to 
soft-inclusive processes in PYTHIA 6 is somewhat crude. 
It uses parametrized cross sections to predict the rates of 
single (SD) and double (DD) diffractive dissociation dif- 
ferentially in the mass(es) of the diffracted system(s) ifTSll . 
Each diffractively excited system is represented by a single 
"string" of the given mass, which is hadronized according 
to the Lund string fragmentation model lfT9ll20ll . We note 
that this type of diffractive modeling can be characterized as 
"soft" since it does not include a mechanism for hard, high- 
mass diffraction, such as diffractive jet production. We in- 
clude it, nonetheless, to give an idea of how the bulk of soft 
diffractive processes affect our conclusions. We also note 
that "typical" minimum-bias cuts are designed to reduce the 
contamination by diffractive processes, such that our con- 
clusions should not depend too crucially on the modeling 
of this component. 

The modeling of inelastic, non-diffractive processes is 
more sophisticated and is based on a picture of multiple 
parton interactions (MPI). In PYTHIA 6, there are two ba- 
sic MPI frameworks available, which we shall refer to as 
"old" [21 1 and "new" ||22l|23l. (The latter is similar to the 
modeling in PYTHIA 8.) Briefly stated, the main differences 
between the old and new models are: 

- Old: virtuality-ordered parton showers, no showers off 
the additional MPI, and a relatively simple description 
of the fragmentation of the beam remnants in which the 
baryon number is carried by the remnant. 

- New: transverse-momentum-ordered parton showers, in- 
cluding showers off the additional MPI, and a more so- 
phisticated treatment of the beam remnant, in which "string 
junctions" 12A\ carry the beam baryon number. 

In both cases, the fundamental MPI cross sections are de- 
rived from a Sudakov-like unitarization/resummation of per- 
turbative QCD 2 — ?• 2 scattering ||2T1 . normalized to the 
total inelastic non-diffrative cross section, and regulated at 
low by a smooth dampening factor. The latter is inter- 
preted as being due to colour screening, and the dampen- 
ing scale, p^Q, represents the main tunable parameter of the 
model. Two other significant parameters are the assumed 
transverse shape of the proton (lumpy or smooth), and the 
strength of colour reconnections (CR) in the final state, cf. lU. 

2.1 PYTHIA Tunes 

We shall consider a small selection of recent tunes that use 
the "old" and "new" models, as follows. Field's Tune DW 
ll25l is currently the "preferred" Tevatron tune and has been 
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extensively tested there. Perugia O also attempts to give a 
good fit to Tevatron data, but uses the new model and incor- 
porates an updated set of fragmentation parameters tuned to 
LEP data by the Professor collaboration ||26l. That col- 
laboration's own tunes of the old and new model are called 
Q20 and PTO, respectively ll26ll . Finally, we also include a 
tune called ACR ||27]| . which represents a hybrid between 
the two models; it uses the basic shower and MPI frame- 
work from the old model with the colour-reconnection (CR) 
model of the new one. It is included to make it possible 
to isolate whether specific features are due only to the CR 
model or not. All tunes were run with Pythia version 6.4.21 
Table [1] show the tunes used along with their three-digit 
codes in the Pythia subroutine PYTUNeB 



Parameter 


DW 


ACR 


Q20 


PO 


PTO 


PYTUNE 


103 


107 


129 


320 


329 



Table 1: Pythia tunes and corresponding three-digit PYTUNE 
codes. 



2.2 Sub-Process Samples 

Pythia includes four distinct processes in its simulation 
of soft-inclusive physics: elastic scattering, single diffrac- 
tive dissocation, double diffractive dissociation, and inelas- 
tic non-diffractive (low-p^) interactions. The sum of these 
contributions is the total hadron-hadron cross section. 

Elastic scattering occurs when the colliding protons in- 
teract without either of the beam hadrons breaking up. We 
shall not consider this source further, since it does not pro- 
duce any particles at central rapidities. 

In single diffractive dissociation (SD), one of the incom- 
ing hadrons breaks up, and the other does not. In this sit- 
uation, a spread of low-px particles is expected from the 
disintegrated system over a limited rapidity region near the 
dissociated hadron, while the undissociated one continues 
with a modified momentum. 

Double diffractive dissociation (DD) involves the break- 
up of both beam particles. Here, both systems generate sig- 
nificant low-p^ particle deposits from disintegration over 
rapidity, typically with a gap between them. 

Low-p^ or non-diffractive interaction involves partonic 
scattering processes, all the way from soft to hard, with 
the latter mapping smoothly onto the dijet tail. In this case 
particle production is more localized, with higher-px con- 
stituents and the possibility — switched on by default — of 
additional perturbative activity such as parton showers and 
multiple parton interactions. 

For all tunes, we start from an inclusive sample com- 
posed of the three inelastic process types, distributed ac- 
cording to their relative cross sections, which are fixed by 

' Note: these tunes can also be activated by setting the parameter 
MSTP ( 5 ) to the relevant PYTUNE code. 



Pythia's default parametrizations ||28]| . Since the descrip- 
tion of the diffractive components is quite simple, it would 
not make much sense to attempt to isolate individual con- 
tributions to the particle production within the two diffrac- 
tive samples. The particle production in the low-p^ sample, 
however, receives contributions from several different algo- 
rithmic components which we may label as: "hard" scatter- 
ing, parton showers, MPI, and remnant fragmentation, each 
with its own distinct behaviour. 

In order to isolate what happens as each component is 
"switched on", we consider four different variants of the 
low-px sample: 

1 . loy/-p± : Everything on, corresponding to the most phys- 
ical description of the low-p^ sample. 

2. HARD: No parton showers, no MPI. I.e., a single par- 
tonic 2—7-2 interaction, with no parton showers. 

3. RAD: No MPI. I.e., a single partonic 2 — )• 2 interaction, 
with parton showers. 

4. MPI: No parton showers. I.e, multiple parton-parton in- 
teractions, without showers. 

Note that all variants are passed through the string fragmen- 
tation model in order to give final-state hadrons. 

For each sample (and for each variation of the \ow-p± 
one), 100,000 collisions were generated at = 900 GeV. 
This sample size is sufficient to overcome statistical fluctu- 
ations for the measurements of interest, and more than the 
data set in |[29l . Obviously, larger sample sizes would still 
be of interest, to study the tails of distributions, but we shall 
here mostly be concerned with the bulk of the physics. 

3 Selection Procedure 

Our starting point is a very loose selection, requiring at least 
some activity in a rather broad central region consistent with 
the charged-particle acceptance region of the UA5 experi- 
ment. Thus, only stable charged particles within a pseudo- 
rapidity range of r/ = (—5., 5.) are selected. Although the 
central trackers of the ATLAS and CMS experiments only 
extend to pseudorapidities of ±2.5, we note that the For- 
ward Multiplicity Detector (FMD) in the ALICE ex- 
perimented could be used to extend charged-particle mea- 
surements of the types we consider to include pseudorapidi- 
ties up to r/ 5. By stable charged particles, we mean all 
charged particles with proper lifetimes cr > 10 mm (thus, 
e.g., A and K hadrons are stable). By default, we do not 
apply any cuts in p^ unless explicitly stated otherwise. Se- 
lected events must have at least one charged particle in the 
?7-range. 

Table |2] shows the percentages of selected events with 
this mild requirement. For elastic events, the corresponding 
percentage would be zero, as the scattered protons continue 
on "down the beam-pipe", outside the range of selection. 
Of the included sub-processes, a significant fraction of SD 
events are rejected already at this point, as are about half 
that fraction of the events labeled DD. These correspond 
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mix 


SD 


DD 


low-pf 


HARD 


RAD 


MPI 


DW 


73% 


70% 


86% 


100% 


100% 


100% 


100% 


ACR 


73% 


70% 


86% 


100% 


100% 


100% 


100% 


Q20 


73% 


69% 


86% 


100% 


100% 


100% 


100% 


PO 


73% 


68% 


85% 


100% 


100% 


100% 


100% 


PTO 


73% 


68% 


85% 


100% 


100% 


100% 


100% 



Table 2: Selection efficiencies for various Pythia sub-processes 
and particle production mechanisms for each tune. Each sam- 
ple consisted of 100,000 generated events. 



to events where the fragments of the dissociated proton(s) 
were concentrated at high rapidities, beyond the acceptance 
region. The events in the DD sample are less likely to be re- 
jected, since they have two "chances" to produce particles 
in the central region. The low-n^ ones are selected with 
approximately 100% efficiencjo. The "mix" sample con- 
tains the sum of all three subprocesses, weighted according 
to their respective cross sections as given by Pythia. The 
particle production samples (hard process only, radiative, 
MPI and combined i.e. low-p^) have similar selection rates 
as they all include a central 'hard' interaction. 

An illustration of the effect of the event selection on the 
inclusive p±^ distribution is shown in figure [T] for the low- 

sample. The top pane shows the p±^ distribution of all 
generated tracks in black and of the selected ones in green, 
for the DW tune. The main effect is a reduction in the total 
number of accepted tracks by 10-15%. A secondary effect, 
however, is a model-dependent hardening of the spectrum. 
We highlight this in the lower pane, which shows the ratio 
of selected to unselected tracks for each of the tunes. We 
observe that the DW and Q20 tunes of the "old" model ex- 
hibit an approximately constant value of this ratio, indicat- 
ing that the shape of the pj^ spectrum is not greatly different 
at high rapidities than in the central region. By contrast, the 
tunes of the "new" model (PO and PTO) exhibit a noticeable 
shape, indicating that for those models, the spectrum of the 
unselected high-rapidity tracks is systematically softer than 
in the central region. Also, since the hybrid ACR tune fol- 
lows the "new" tunes here, we may interpret this behaviour 
as being related to the CR model used in the new tunes. 

We note that these differences would translate directly 
to uncertainties for any purely model-based extrapolation 
beyond the central region. For the remainder of this report, 
however, we shall be concerned only with the central region 
itself. 



4 Inclusive Distributions 

Current LHC studies of minimum-bias events, e.g. ||9]- [T7]| . 
have focused mainly on the "basic four" charged-particle 
distributions: P{nch), dnch/drj, drich/dpt, and (pt) vs. rich- 



For completeness, we note that a few of the generated low-px 
events do fail, below the per mille level, having produced no or only 
neutral particles in the central region. 
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Figure 1: Above: Charged particle p± distribution for low-p^ 
sample of DW tune, with (green) and without (black) selec- 
tion cut. Below: Ratio of selected to unselected tracks for the 
low-p± sample of each tune. 



In this section, we comment briefly on these and illustrate 
the behaviour of our chosen set of Monte Carlo tunes for 
later reference. More comments on these distributions can 
be found, e.g., in Ed. 



4.1 Multiplicity 

In table [3j we compare the average charged particle multi- 
plicity, {rich), between tunes, for the three different inelastic 
samples, as well as for their cross-section-based mixture. At 
this point, we subject the samples only to the selection re- 
quirement \rj\ < 5 mentioned above. Table |4] contains an 
equivalent comparison, for the different physics variations 
of the low-p^ sample. To help illustrate the overall spread 
in predictions, we also quote a "range" of variation, at the 
bottom of each table, which is defined as the highest aver- 
age multiplicity of the tunes minus the lowest, normalized 
to the lowest multiplicity, i.e. 



range = {{rich) 



{nch)^J/ {nch). 



(1) 



The range of {rich) predicted within jr/j < 5 varies by 
10-20% between the different tunes for each of the inelas- 
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Tune 


SD 


DD 


low-px 


mix 


DW 


10.1 


12.0 


38.6 


23.5 


ACR 


10.1 


12.0 


38.4 


23.5 


Q20 


10.3 


12.1 


38.3 


23.4 


PO 


9.1 


10.8 


35.4 


21.7 


PTO 


9.1 


10.8 


36.8 


22.5 


range in (tic/j) (%) 


14.1 


12.8 


9.0 


8.5 



Table 3: Average charged particle multiplicity for the various 
minimum-bias sub-processes and tunes. 



Tune 


HARD 


RAD 


MPI 


All on 


DW 


28.0 


32.9 


33.5 


38.6 


ACR 




31.6 


36.2 


38.4 


Q20 


29.1 


30.9 


36.3 


38.3 


PO 


23.8 


26.0 


29.6 


35.4 


PTO 


24.1 


26.0 


31.1 


36.8 


range in (nch>(%) 


22.4 


26.6 


22.4 


9.0 



Table 4: Average charged particle number for each particle 
production mechanism. Note: the column labeled "All on" is 
identical to the "low-px" one in table |3] 

tic sub-processes in table [3] For the diffractive processes, 
there is no parton showering and no MPI. The considerable 
differences between models are therefore solely generated 
by the different tunings of the hadronization model. Since in 
particular Q20 and PTO were tuned to exactly the same LEP 
data by exactly the same tuning program (PROFESSOR), the 
difference between them here highlights the need for in situ 
constraints on the non-perturbative fragmentation parame- 
ters. It also indicates the need for a Monte Carlo modeling 
of diffraction that would be more theoretically consistent 
with the treatment at LEP, where the fragmentation param- 
eters are dependent, e.g. on the perturbative parton shower 
cutoff. For the time being, we conclude that the fragmenta- 
tion tuning of the new model used by PTO and PO produces 
fewer particles in and of itself than that of the old mod- 
els. This appears to be consistent with the distributions 
shown in[3j fewer particles produced from the same energy 
of collision will result in a larger proportion of high-p^ par- 
ticles in selected events. Though not the main focus of this 
report, this would clearly be worth a more detailed analysis 
especially in the context of diffractive studies. 

4.2 Track p± 

Figure |2] shows the distributions for the low-p^ sample 
of each tune, with the range of variation of the average {p±^) 
spanned by the tunes for all the sub-samples given in table 
|5j with the range defined as in eq. ([1]). 

As before, we note that the tunes exhibit differences of 
the order 10-20%. In the low region (p^ < 1 GeV) the 
old shower tunes DW, ACR and Q20 lie above the new 
shower models PO and PTO. This reverses in the region 
above 2 GeV. Due to limited statistics, we do not plot the 
tail of very high-p^ charged particles here, but note that the 



10-'^ 




P, [MeV] 

Figure 2: Selected charged particle logarithmic p± distribution 
for low-pj^ sub-sample Pythia tunes. Lower pane: ratio to the 
DW distribution. 





mix 


SD 


DD 


low-pf 


HARD 


RAD 


MPI 


range in (pt) (%) 


5.7 


3.1 


3.5 


6.1 


10.5 


6.1 


10.5 



Table 5: Range of mean pt values for various Pythia sub- 
processes and particle production mechanisms for each tune. 

trend of the new models to generate harder p±_ tails is illus- 
trated in 111. 

4.3 Track rj 

Figure [3] shows the 7] distributions for the selected mod- 
els. Differences are again at the 10-20% level. Especially in 
the context of this distribution, such differences have some- 
times been represented as "large". Let us recall, however, 
that the PYTHIA modeling is rooted in perturbative QCD, 
that we are here dealing with processes which have no hard 
scale, and that the number of charged particles is not an 




Figure 3: Selected charged particle logarithmic 77 distribution 
for low-pj^ sub-sample PYTHIA tunes. Lower pane: ratio to the 
DW distribution. 
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infrared safe observable. In this light, while we might still 
hope to constrain the modeling better, we nonetheless wish 
to point out that it is, in our opinion, grossly misleading to 
characterize order 10% differences as large. 

Indeed, the small differences between tunes are high- 
lighted by the zero-suppressed Y-axis in the plot. Thus, while 
there is clearly some sensitivity to central vs. forward pro- 
duction mechanisms in this distribution, its ability to dis- 
criminate between models is still limited. Agreement be- 
tween each tune is generally good, especially in the most 
easily observable region, \r]\ < 2.5. We conclude that addi- 
tional, linearly independent, information on the structure of 
events in rj, could provide valuable additional constraints. 

5 Forward-Backward Correlations 

We come now to the main part of this report, in which we 
study several types of forward-backward correlations, b, for 
different production mechanisms, cuts, and correlation re- 
gions. 

The purpose of these distributions is to enhance the dis- 
criminating power between models, and to reveal their prop- 
erties more clearly, as compared to what can be achieved 
with the list of observables discussed in section |4l In par- 
ticular, the collinear singularity structure of bremsstrahlung 
corrections in perturbative QCD causes initial- and final- 
state shower activity to generate strong but primarily short- 
range correlations, spanning at most a few rapidity units, 
whereas coloured exchanges between the beam hadrons (e.g., 
MPI) can generate correlations that are weaker but which 
span the entire rapidity range between the remnants. Thus, 
the shapes and normalizations of the b distributions contain 
valuable information on the relative dominance of differ- 
ent particle production mechanisms, information which we 
argue is linearly independent from that contained in the cur- 
rent "standard" distributions. 

This section is divided as follows: we first consider a 
standard inclusive "minimum-bias" b correlation in section 
15.11 illustrating how it is affected by different choices of bin 
size and by pj^ cuts; in section [5^ we illustrate the sensitiv- 
ity of this correlation to different particle production mech- 
anisms, using the HARD, RAD, and MPI samples defined 
above, and to contamination by diffractive processes (SD, 
DD). In this way, we gain a map of how different cuts and 
different process mixtures affect the correlations, that we 
hope will be useful for future reference. We shall seek to 
extract further information by defining also a set of b cor- 
relations that are sensitive to the azimuthal structure of the 
events, which will be the focus of section 15.31 We shall re- 
fer to these latter observables, which are essentially binned 
double -differential r]-(j) correlations, as "twisted" b correla- 
tions. 

5.1 Inclusive b Correlation 

The standard b correlation is defined as: 




n 

Figure 4: b correlation for selected events with various pseudo- 
rapidity bin sizes, Arj. The single-point correlation for a bin 
size of Z\?/=5.0 (red) is shown covering the whole ?/-region. 



(T(nfc)cr(n/) - {uff 

where Uf (rih) is the activity in a specific forward (back- 
ward) region of the detector. "Activity" can be measured 
by a number of observables in the detector, e.g., energy, 
charged particle multiplicity (inclusively or above a given 
Pj^ threshold), momentum sum, etc. Here, we shall focus 
on the charged-particle multiplicity, as has also been done 
in most previous studies, though we emphasize that, e.g., 
calorimetric energy sums, could also be interesting to ex- 
plore (see, e.g., Q). 

The "forward" and "backward" regions are defined by 
bins of a specific size in — typically chosen to be be- 
tween 0. 1 and 1 .0 unit wide inrj — which are separated by 
some variable distance and arranged symmetrically around 
a midpoint which is usually taken to be the centre of the 
detector, rjc = 0, corresponding to the CM of the colliding 
hadrons. Although we shall not do so here, we note that cor- 
relations between the central and forward region are also of 
interest and can be probed, for example, by fixing one bin in 
the central region and letting the other slide into the forward 
or backward region, corresponding to choosing rjc ^ 0. A 
study somewhat along these latter lines has been performed 
by UA5 ||29l and could also be interesting to maximize us- 
age of the asymmetric coverage of the ALICE FMD. 

The optimum bin size to use in eqn. ([Jll is a function of 
statistics and of the r/-range observed. If the bin size is too 
small, genuine correlations will be washed out by statistical 
fluctuations. With too large a bin size, the resolving power 
of the correlation over the limited ?7-range will be lost. 

Figure |4] contains a comparison of the b correlation vs. 
r] (specifically the 7] value of the centre of the forward bin, 
with the backward one located at —7]) for varying bin sizes 
from 0.03 to 5 units wide, without imposing any p±^ cuts 
at this point. Obviously, the largest sizes are too coarse to 
discern much structure in the correlation distribution. Mid- 
range bin sizes, Ar]= 1.0,0.5 and 0.25 exhibit best the trends 
over the r/-range; for this particular model (DW), a high cor- 
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Figure 5: The b correlation distributions for hadron-level 
charged particles for the DW tune with various explicit p± cuts. 
Lower pane: ratio to the non-cut distribution. 
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Figure 6: The b correlation distributions for hadron-level 
charged particles for the Perugia tune with various explicit 
p± cuts. Lower pane; ratio to the non-cut distribution. 



relation at low t] can be distinguished from a mid-rj plateau 
and a further drop in correlation at high ij. Going to even 
smaller bin sizes, Ari=0.l25, 0.0625 and 0.03125, we be- 
gin to lose the structure in the distribution as statistical fluc- 
tuations start to dominate. We conclude that a bin size of 
Ai]=0.5 is reasonable for this study. Note, however, that go- 
ing to different CM energies and/or imposing cuts could 
change this conclusion; the average accepted multiplicity at 
each r] value determines the relative size of the statistical 
fluctuations and hence affects the optimum bin size. As a 
consequence, it is not possible, therefore, to directly com- 
pare b distributions taken with different cuts, or which use 
different-sized Ar] bins. 

To illustrate the effect of p±^ cuts, figs. |5] and |6] show the 
b correlations subjected to a range of different cuts for 
a tune of the old (DW) and new (Perugia 0) model, respec- 
tively. (Note that these cuts are applied also at the level of 
the event selection, so only events with at least one particle 





central bin 


mid-range bin 


extreme bin 


Tune 


bo 


6()„5/6(. 


6i,.!/6o 


bi, 


60.5/60 


61.5/60 


61] 


60.5/60 


61,5/60 


DW 


0.38 


0.79 


0.39 


0.28 


0.74 


0.24 


0.06 


0.46 


-0.08 


ACR 


0.45 


0.71 


0.34 


0.28 


0.68 


0.30 


0.05 


0.38 


-0.03 


Q20 


0.43 


0.74 


0.29 


0.27 


0.68 


0.17 


0.05 


0.42 


0.16 


PO 


0.46 


0.72 


0.23 


0.23 


0.66 


0.18 


0.00 


0.71 


-0.39 


PTO 


0.45 


0.73 


0.23 


0.24 


0.65 


0.18 


0.01 


0.15 


-0.05 



Table 6: Central(T7=0-0.5), mid-range (77=2.5-3) and extreme 
(77=4.5-5) correlation values with fraction of correlation re- 
maining after p^-cut= 0.5GeV and 1.5 GeV for various Pythia 
tunes at \/s=900GeV. 60 is the coiTelation value for p±- 
cut=O.OGeV, 60.5 the value with pj^-cut=0.5GeV and 61.5 the 
value p^-cut=1.5GeV. 

harder than the given cut are included, for each curve.) 
One clearly sees the reduction of the correlation strengths 
as fewer tracks make it into each bin. Note that only posi- 
tive correlations are expected and plotted. We interpret any 
small negative correlations as arising from statistical fluc- 
tuations in poorly populated bins. 

Table |6] summarizes the effect of p±^ cuts by giving the 
b correlation values in the central (| 7/1=0-0.5), mid-range 
(|r7|=2.5-3) and extreme (|?7|=4.5-5) bins, without any pj^ 
cut, together with the reduction in the correlation strengths 
caused by pj^ cuts of 500 MeV and 1.5 GeV. In the cen- 
tral r^-region, the effect of the p±^ = 500 MeV cut (pink 
dashed line in figs. |5] and IS is to lower the correlation by 
20 - 30%. When the p± cut is increased to 1.5 GeV (dark 
purple solid line, second from bottom in the plots), the re- 
duction is much more severe, and more interestingly is not 
a simple scaling from the reduction caused by the previ- 
ous cut. Thus, e.g., ACR exhibits the largest reduction with 
the first cut, but the second-to-least with the second. Also 
interestingly, this pattern changes as one goes from mid- 
range to the extreme rj bins. (At the extreme end, though, 
the correlations are small, especially with high-pj^ cuts, and 
hence are easily overpowered by statistical fluctuations.) 
When comparing correlation values surviving the 500 MeV 
cut between central and extreme ?7-bins, we may conclude 
that the particle momentum distributions are heterogeneous 
across the r^-range, due to the different parameters between 
tunes and differing particle production mechanisms. Mea- 
surements should therefore by no means be restricted to the 
most inclusive definition possible for a given experiment. 

As a final summary of this part of the study, fig. |7] con- 
tains a comparison of the b correlation distribution for each 
of the different tunes, again without p±^ cuts imposed. As 
was already apparent from figs. [5] and |6l the old and new 
models exhibit qualitatively different shapes. We interpret 
this in the following way: due to the inclusion of showers 
off the MPI in the new models, more of their total parti- 
cle production is driven by shower activity than what was 
the case in the old ones, which have a larger average num- 
ber of MPI i3j. The new models therefore exhibit stronger 
short-range correlation^ and weaker long-range ones than 

Note that these particular tunes of the new model have fewer aver- 
age charged particles than those of the old, cf. table[3] Due to the dilu- 
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Figure 7: Inclusive b coiTelation distribution for selected low- 
p± Pythia tunes. Lower pane; Ratio to the DW distribution. 



their older counterparts, with a crossover point somewhere 
around Ir^j = 1 — 2. 

There are further differences in the b correlation dis- 
tributions especially within the old model class. DW has 
the most distinctive shape of the old models, with a clear 
plateau-like structure at mid-r] which is not as pronounced 
for any of the other models. This is consistent with the 
dN/drj distribution being higher for this tune for \r]\ > 3 
than for any of the other models, cf. fig.[3j Due to this sig- 
nificant shape difference, the Q20 distribution, for instance, 
lies above DW at low rj, but then drops below it at high i]. 
It should therefore be clear that a measurement of the shape 
of this distribution out to as high r] as possible would yield 
valuable information. 



5.2 Physical Sources of Correlations 

To investigate the sensitivity to the different sources of par- 
ticle production in more detail, we now turn to the HARD, 
RAD, and MPI samples, as compared to the default low- 
sample which has all the physics components switched 
on. The results for one tune of the old model (DW) and 
one of the new (Perugia 0) are shown in figs. [8]and|9j re- 
spectively. In both cases, and also for the other tunes not 
shown here, the general trend is for the MPI component to 
dominate the distributions. Again, this has partly to be un- 
derstood in the light of the MPI component generating the 
largest part of the multiplicity, see table |4l such that sta- 
tistical fluctuations are relatively more important when that 
component is switched off, as in the RAD and HARD sam- 
ples. 

As discussed above, however, one notes that the HARD 
component by itself only produces very short-range corre- 

tion effect caused by statistical fluctuations, their correlation strengths 
are therefore intrinsically a bit lower than if they had been made to 
give the same average multiplicities as their older counterparts. 




Figure 8: Inclusive b correlation distribution for tune DW parti- 
cle production mechanisms: low-p±, hard process (HARD), ra- 
diative production (RAD) and multi-parton interactions (MPI). 
Lower pane: ratio to the low-p± distribution. 
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Figure 9: Inclusive b correlation distribution for tune Peru- 
gia particle production mechanisms: low-p^, hard process 
(HARD), radiative production (RAD) and multi-parton inter- 
actions (MPI). Lower pane: ratio to the low-pj^ distribution. 

lations, that drop off quickly to zero. Adding parton show- 
ers, in the RAD samples, extends the reach of these correla- 
tions somewhat further in 7], including a small tail towards 
very large r], presumably generated by initial-state radiation 
from the beams. 

Interestingly, the behaviour of the MPI component is 
somewhat different between the two kinds of models. In the 
old model, fig. [8l the MPI component becomes completely 
dominant at large r] and there has the same magnitude as 
the low-px sample itself. In the new model, fig. |9l the MPI 
component alone drops off and is eclipsed by the shower 
component at the highest values of r], indicating a qualita- 
tive difference between the models, consistent with the new 
model deriving more of its total particle production from 
shower-related activity. 

The effects of diffractive components, a non-zero con- 
tamination of which may be present especially in very in- 
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Figure 10: Inclusive b correlation distribution for tune DW 
minimum bias sub-processes. Lower pane; ratio to the distri- 
bution of the low-px sample. 

elusive minimum-bias measurements, are illustrated in fig. 
[TOl The correlations in the SD and DD samples are intrinsi- 
cally shorter-range than those of their non-diffractive coun- 
terparts, consistent with diffractive systems having a limited 
extension in rapidity. 

However, we also see an interesting effect of combining 
the samples, namely that the correlations in the combined 
sample are stronger than in any of the individual compo- 
nents. It seems that, by mixing in less correlated diffrac- 
tive components, we have actually enhanced the final cor- 
relations. What is going on? This effect can be illustrated 
by imagining we have two separate distributions, A and B 
(in our case represented, e.g., by the diffractive and non- 
diffractive samples). Imagine further that the fluctuations 
inside each sample are purely statistical, for illustration, 
such that the correlation strength inside each sample is zero. 
What will happen when we look at the combination A-\-Bl 
If the mean of A is smaller than that of B, then every event 
of type A will look like it fluctuated down, systematically, 
from the mean of A + B, and conversely for the B sample. 
In their combination, therefore, we will see a non-zero cor- 
relation if the mean values are different. Since the diffrac- 
tive and non-diffractive event samples have very different 
average multiplicities (see table [3]l this effect will lead to 
an increase in the correlations in the combined sample, as 
observed in fig. \T0\ 

5.3 'Twisted' b Correlations 

We now turn to the dependence on azimuth of the forward- 
backward correlation strengths. A related type of correla- 
tions sensitive to both rj and azimuthal (j) were recently 
highlighted by the CMS experiment ||3T1 and have stimu- 
lated quite a lot of interest due to the observation of the so- 
called "ridge effect" in high-multiplicity events. The cor- 
relations presented in this report are somewhat simpler in 
spirit, and our focus is not primarily on high multiplicities. 



but we note that it could be an interesting follow-up study 
to determine whether twisted 6-correlations could also be 
used to shed more light on the ridge. 

We shall study the (p dependence of the 6-correlations in 
two ways. The first is based only on the detector geometry. 
As this is independent of the event shape, no preference is 
given to any particular direction. The second method gives 
preference to the direction of the leading charged particle 
in the event. This will bias the zero point in (j) to coincide 
with the most active part of the event. 

In each case, we divide the c^-plane into three regions 
of size A(j) = 2 /Sit. For the detector-defined geometry, 
we define a parallel, an opposite, and a transverse region. 
(Note: we use "parallel" and "opposite" here, to distinguish 
the geometry from Field's "towards" and "away" regions, 
the latter of which we take to be defined relative to the di- 
rection of a lead particle or jet and not by the absolute de- 
tector geometry.) Quite arbitrarily, we define 

- A/7ara//e/ region covering the region —vr < (f) < — l/Svr 
in absolute azimuth, 

- An opposite region covering < (p < 2/3tt, 

- A transverse region occupying the region between these, 
i.e., the slices — l/37r<0<O and 2/3tt < (j) < tt. 

In calculating the correlation between ij — (j) regions the 
comparison is always to the parallel case on one side. We 
are aware that this is quite crude and that one could increase 
statistics by integrating over the location of the arbitrary az- 
imuthal zero point, but point out that this is intended merely 
as a first exploration of the properties of 'twisted' correla- 
tions. 

The terms of the b correlation expression now refer to 
rj-hins with a i;^>-dependence. Hence, the correlation expres- 
sion must include this new degree of freedom. Since all re- 
gions are a priori equivalent, the normalizing terms in b, 
{uf)" and (rij:), are taken simply from the parallel one. 
Only the product of activity in corresponding bins of r] — (p 
are sensitive to the variation in cp region. The new expres- 
sion, 6*^^"^', for the correlation becomes: 




Figures [TT] and [12] show the three different types of cor- 
relations we can obtain for the detector-based geometry, 
for Pythia tunes DW and Perugia 0, respectively. As ex- 
pected, both models exhibit a peak in the correlation at low 
r] in the parallel region, illustrating that the low-?] corre- 
lation is also most pronounced at low A<p. Beyond the first 
bin in 7], however, the opposite correlation is strongest. This 
follows from momentum conservation, smeared out over rj. 
The correlation with the transverse region is as close as we 
can come to defining an "underlying event" in an otherwise 
featureless minimum-bias event without a reference direc- 
tion. 
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Figure 11: 6^™'^* correlation distributions for hadron-level 
charged particles for DW Pythia tune for the three different 
combinations of regions, defined with respect to the absolute 
detector geometry. Lower pane: ratio to the parallel distribu- 
tion. 
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Figure 13: correlation distributions for hadron-level 

charged particles for the DW tune for the three different com- 
binations of (p regions, defined with respect to the lead particle 
trajectory. Lower pane: ratio to the towards distribution. 
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Figure 12: fo^^'^* correlation distributions for hadron-level 
charged particles for the Perugia tune for the three different 
combinations of (p regions, defined with respect to the absolute 
detector geometry. Lower pane: ratio to the parallel distribu- 
tion. 

The difference in correlation strength between the three 
regions are not extremely large in absolute terms, however. 
This leads us to consider whether there is a way to en- 
hance the differences while remaining in a minimum-bias 
context. By choosing the zero point of the (j) coordinate, 
event by event, to be the direction of the leading (hardest) 
charged particle, we can now use the nomenclature of Field 
and define the towards region to include the azimuthal an- 
gles (j) < zbl/Svr around the lead particle, the away region 
covers ±2 /Sir < (j) < ztyr and the transverse region lies 
in between, over ibl/37r < (j) < ib2/37r. We label the 
corresponding correlation b^^^'^, to distinguish it from the 
detector-based one. 

This orientation has a significant effect in particular for 
events with a semi-hard perturbative scattering, in which 




n 

Figure 14: 6'^^'^ correlation distributions for hadron-level 
charged particles for the Perugia tune for the three different 
combinations of cj) regions, defined with respect to the lead par- 
ticle trajectory. Lower pane: ratio to the towards distribution. 



the main axis in the transverse plane becomes oriented to 
the production axis of the outgoing partons. The bias to- 
wards (/> = as the direction of the lead particle means that 
the three different (p regions can no longer be expected to 
have the same averages and variances. Nonetheless, in or- 
der to define a measure comparable to the one above, we 
shall still define the normalizing terms with respect to the 
towards region, so that eq. ([3]) still holds, although its sta- 
tistical interpretation is modified. 

Redoing the twisted correlations with this definition of 
the zero point in azimuth, we obtain figures [13] and [141 for 
DW and Perugia 0, respectively, with the other tunes ex- 
hibiting similar qualitative features. The general remarks 
are similar to those for the detector-based geometry, but the 
differences between the regions are much more clearly vis- 
ible. Also, the transverse region can here be clearly iden- 
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tified as lower than the others, consistent with its being an 
"underlying event" to the production of a "hard particle". 



6 Event Shapes 

A further characteristic of the structure of the events in the 
minimum-bias sample can be gained by considering their 
eigenvalues along the principal event axes in the transverse 
plane. Obviously, we expect minimum-bias events to be 
much more uniform than, e.g., jet events, and hence to have 
smaller eigenvalues, but how much smaller? What is their 
average shape, and how much does it fluctuate between 
events? 

To address these questions, and to gain a first idea of 
their sensitivity to the physics modeling, we consider the 
transverse thrust (T^) and transverse minor (Mt) values and 
axes 1321 . 
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Figure 15: Transverse thrust distributions for the low-px sub- 
sample of selected tunes, normalized to unit area. Lower pane: 
ratio to the DW distribution. 



6.1 Transverse Thrust 

The transverse thrust axis can be found by maximizing the 
coincidence of an arbitrary vector with the dominant direc- 
tion of particle flow in an event in (j). The value of transverse 
thrust is then defined as: 



max 



(4) 



interesting that, in the old model (DW), the MPI component 
by itself (solid blue lines with triangular symbols) only re- 
duces the peak value very slightly, whereas the addition of 
radiation (RAD: solid green lines with square symbols) pro- 
duces a much larger shift. In the new model, however, the 
MPI and RAD samples each appear to give a similar-size 
shift. Despite their apparent similarities, there are therefore 
still interesting differences underlying these distributions, 
which, as we have argued, the measurement of b correla- 
tions can help resolve. 



where i runs over the charged tracks in the event, is 
the transverse thrust axis unit vector and pY the track 
transverse momentum vector. The observable is bounded 
by 0.5 < < 1.0. Dijet-like events, where the highest- 
momentum particles are produced back to back, have a pen- 
cil-like shape in 0, with particle production aligned pre- 
dominantly along the axis. Such events have high transverse 
thrust values ~ 1. In contrast, in events where non-pertur- 
bative and/or MPI production is predominant, more parti- 
cles will lie off the main production axis, giving a more 
circular distribution of tracks in (p, for which the transverse 
thrust value will lie closer to 0.5. 

Fig(T5] shows the transverse thrust distributions of the 
\ow-p± sub-samples of the selected tunes. Somewhat sur- 
prisingly, perhaps, the models agree to within 10-20% over 
most of the range. This presumably reflects the fundamen- 
tal similarity between the MPI-based perturbative modeling 
in these tunes. A comparison with experimental data could 
give valuable insights on whether the real world is more or 
less "jetty" than these model predictions indicate. 

To further analyze the structure of this distribution, figfT6l 
shows the transverse thrust distributions for the HARD, RAD, 
and MPI samples, compared to the full low-pj^ simulation, 
for the DW (above) and Perugia (below) tunes. For the 
HARD sample (i.e., before showering and MPI), the distri- 
butions are more pencil-like, peaking at somewhat higher 
values of Tj_, illustrated by the dashed (cyan) curves. It is 



6.2 Transverse Minor 

The transverse minor axis lies perpendicular to the thrust 
axis in 6. It is defined as: 



Ml 



max 



(5) 



using the same definitions as in eq. FigUS] shows the 
transverse minor distributions of the low-p^ sub-samples of 
the selected tunes. As before, the variations between models 
is relatively mild, with the contributions from each model 
component, fig. [181 exhibiting similar differences as for T^. 



7 Conclusion 

We have illustrated that forward-backward correlations can 
be used to extract information on the relative strengths of 
different sources of particle production in minimum-bias 
events: models dominated by a single hard (dijet) interac- 
tion exhibit strong short-range correlations and weak long- 
range ones, while models with a larger component of soft 
production between the remnants generate stronger long- 
range correlations. 
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Figure 16: Transverse thrust distributions for the low-pj^ sam- 
ple for a) DW and b) Perugia 0, normaHzed to unit area. Lower 
panes: ratio to the respective low-p± samples. 



We propose to add these distributions to the "standard" 
ones already measured by the LHC experiments, and fur- 
ther to add correlations between different regions in az- 
imuthal (j), which we label 'twisted' forward-backward cor- 
relations. 

We have illustrated these inferences by comparing a small 
set of recent tunes of the Pythia 6 Monte Carlo model. 
Although they are all based on a picture of multiple par- 
ton interactions (MPI) interfaced to the Lund string frag- 
mentation model, they differ qualitatively in the shower and 
remnant modeling, and quantitatively in the fragmentation 
tuning and amount of showering vs. MPL 

We further believe that measurements of event shapes, 
such as transverse thrust and transverse minor, can help 
shed light on the overall properties and structure of minimum- 
bias events. For instance, a model with a strong dominance 
of perturbative (mini-)jet production should also predict event 
shapes closer to equivalent pQCD ones in dijet events, while 
models characterized by other particle production mech- 
anisms should exhibit spectra further from the factorized 
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Figure 17: Transverse minor distributions for lov/-p± sub- 
samples of selected tunes, normalized to unit area. Lower pane: 
ratio to the DW distributions. 

pQCD prediction. In that context, however, the models stud- 
ied here appear to give relatively similar results, presum- 
ably owing to the significant properties they share at the 
level of the underlying perturbative modeling. 
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